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The tidal disruption of a solar mass star around a supermas- 
sive black hole has been extensively studied analytically 1,2 and 
numerically 3 . In these events the star develops into an elongated 
banana-shaped structure. After completing an eccentric orbit, the 
bound debris fall onto the black hole, forming an accretion disk 
and emitting radiation 4 5 6 . The same process may occur on plan- 
etary scales, if a minor body orbits too close to its star. In the Solar 
System, comets fall directly onto our Sun 7 or onto planets 8 . If the 
star is a compact object, the minor body can become tidally dis- 
rupted. Indeed, one of the first mechanisms invoked to produce 
strong 7 ray emission involved accretion of comets onto neutron 
stars in our Galaxy 9 . Here we report that the peculiarities of the 
'Christmas' 7 ray burst (GRB 101225A 10 ) can be explained by 
a tidal disruption event of a minor body around a Galactic iso- 
lated neutron star. This result would indicate either that minor 



bodies can be captured by compact stellar remnants more fre- 
quently than it occurs in our Solar System or that minor body 
formation is relatively easy around millisecond radio pulsars. A 
peculiar Supernova associated to a GRB may provide an alterna- 
tive explanation 11 . 

GRB 101225A image-triggered the Burst Alert Telescope (BAT) 
onboard the NASA's Swift satellite on 25.776 December 2010 UT. The 
event was extremely long, with a T90 (the time interval during which 
90% of the flux was emitted) > 1.7 ks, and smooth 10 by comparison 
with typical 7— ray bursts 12 (GRBs). The total 15-150 keV fluence 
recorded by the BAT is £ 3 x 10 -6 erg cm -2 and there are no signs 
of decay. The X-ray Telescope and the Ultra Violet Optical Telescope 
on board Swift found a bright, long-lasting counterpart to the GRB. 
Strong variability is observed in the early X-ray light curve. The op- 
tical counterpart, which was detected in all of the Ultra Violet Optical 
Telescope filters, lags the X-ray light curve (Fig. 1). The X-ray and 
optical light curves are reminiscent a shock break-out event observed 
in association with GRB060218 1,i , but are fainter (~ 3.5 mag, i.e. a 
factor of ~ 25) and do not have an X-ray afterglow at later times or 
a bright supernova component. Measurements by the European Space 
Agency's XMM-Newton space observatory failed to detect an after- 
glow with an upper flux limit of ~ 10~ 14 erg cm -2 s _1 at the 3er 
confidence level (0.5-10 keV, observation made at At = 23 d after 
the trigger). 

Ground-based telescopes also followed the event, mainly in the R 
and I bands (Fig. 1). Surprisingly, optical spectra recorded in the first 
few days after the event failed to detect any spectral feature, instead re- 
vealing a smooth blue continuum 14 ' 11 . At later times, the optical light 
curves revealed a color change from blue to red. The NASA Hubble 
Space Telescope imaged the field at At = 20 d, finding a quite red ob- 
ject (Veo6w(AB) — 24.6 and V435W — Bqo6w = 1.6), with no host 
galaxy apparent at magnitudes at most 1 mag fainter 15 . Later observa- 
tions carried out at the INAF Telescopio Nazionale Galileo recorded 
an /—band magnitude of I = 23.64± 0.24 at At = 44 d but recorded 
a non-detection (J > 24.5, 3 a confidence level) at At — 55 d. Obser- 
vations at the Gran Telescope Canarias (At = 180 d) detected GRB 
101225AatgAB = 27.21 ± 0.27 and r' AB = 26.90 ± 0.14 (ref. 11). 

These observational characteristics make GRB 101225A unique. 
Without a-priori knowledge of the distance, this event can result from 
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different progenitors. We explored several scenarios (Supplementary 
Information). Extragalactic models (involving GRBs, supernovae or 
blazars) are unlikely on the basis of total energy arguments, in com- 
bination with the faintness of the optical counterpart and putative host 
galaxy. Galactic models (involving X-ray binaries, magnetars or flare 
stars) are ruled out as well, on the basis of the extremely faint qui- 
escent counterpart and on the relatively high Galactic latitude. The 
timescale of the event, the longer duration at longer wavelengths and 
the faintness of a quiescent counterpart at all wavelengths all suggest 
a disruptive event involving a compact object. We propose that GRB 
101225 A results from the tidal disruption of a minor body around an 
isolated neutron star. Similar models have been already put forward 
to explain the presence of metal absorption lines in the white dwarf 
optical spectra 16 and the infrared excess in the spectra of hot white 
dwarfs 17 . 

A minor body gets disrupted if it comes within a distance of ~ 
10 5 — 10 6 km from the neutron star, where the internal forces that hold 
the body together are overwhelmed by the tidal pull of the neutron star. 
The debris that remain bound are thrown into highly eccentric orbits 
and then fall back to form an accretion disk around the star 18,19 . If the 
minor body has a low tensile strength it gets fully disrupted, its early 
behaviour can be chaotic owing to the formation of the disk with flares 
coincident with the periastron passages 6,20 , and the bolometric light 
curve is expected to decay, in the long term, as i~ 5//3 (ref. 1-2), 

If the matter that falls back to form the disk accretes at a sub- 
Eddington rate at all times (as in our case), the disk emission can 
be described as a spatial superposition of black bodies of increasing 
temperatures and depends on the parameters of the encounter 21 . The 
physics of the event depends, to a first approximation, only on a very 
few parameters: the periastron of the minor body's orbit (r p ), the mi- 
nor body's mass and radius (M* and 7?*; the mean density of a minor 
body in the Solar System is in the range ~ 1 — 10 g cm~ 3 , which 
adds an additional constraint), the minor body's tensile strength, the 
neutron star mass (Mns) and the source distance. 

We fit the X-ray (1 keV), ultraviolet (UVW2 and UVW1 bands, 
centred at 2,030 and 2,634 A, respectively) and optical (R and / 
bands, centred at 6,400 and 7,700 A, respectively) light curves of GRB 
101225 A with a four-parameters phenomenological model of tidal dis- 
ruption onto a neutron star 21 (Fig. 1). We assume that the neutron star 
mass is 1.4 Mq and that the minor body has a density of 1 g cm and 
a low internal tensile strength. Our model is indicative of the physics 
involved, even though it cannot capture the full details of the event. 
We find an adequate fit for M»~ 5 x 10 20 g (half the mass of the 
asteroid Ceres) and a periastron radius of ~ 9, 000 km (which is well 
within the tidal radius, thus demonstrating the internal consistency of 
the model). The mass of the minor body can be estimated within a 
factor of < 3. This mass is large by comparison with those of objects 
in the asteroid belt but is typical of objects in the Kuiper belt. 

Assuming that half of the mass is accreted onto the neutron star 
and the total fluence is 4x 10~ 5 ergcm~ 2 , estimated within the 1 — 10 5 
s time interval (including a bolometric correction motivated by spec- 
tral fitting), we derive a distance of ~ 3 kpc, placing GRB 101225A in 
the Perseus arm of the Galaxy. The observed radiation is emitted both 
by an accretion disk (with an outer radius twice as large as periastron) 
and during the impact of matter onto the neutron star surface. The 
early spectral energy distribution (SED) at five different epochs can be 
adequately fitted by the emission from an accretion disk (with a tem- 
perature of T d ~ 1.8 keV in units of energy, in agreement with model 
predictions) plus the emission from a boundary layer in which the disk 
matter slows down to accrete onto the neutron star (Supplementary In- 
formation). In particular, the initial high-energy emission observed by 



the BAT can be modelled as a black body with T ~ 10 keV 11 . This 
cannot come from the disk but can be easily accounted for by mass 
accretion onto the neutron star surface. The presence of the boundary 
layer component rules out the compact object's being a black hole. 
Small differences in the return times of debris are expected during the 
first periastron passages, and these give rise to X-ray variability on the 
timescale of this fall back 6 ' 20,22 , as observed in the X-ray band (Fig. 
1, inset). The X-ray variability timescale is indeed very similar to the 
fall-back timescale obtained from our modeling. 

The phenomenological model can naturally explain the (observed) 
slower decay at longer wavelengths. It somewhat underestimates the 
observed decay of the optical flux starting at ~ 20 d. Changes in the 
properties of the accretion disk are expected on the basis of the the- 
ory of accretion onto neutron stars. The inner disk radius is increased 
by the neutron star's magnetic field (even when this takes its minimal 
value, 10 8 — 10 9 G). Furthermore, when the accretion rate decreases 
below a critical rate, the disk develops a well-known thermal- viscous 
instability leading to a fast cooling 23,24 . We have calculated that, on 
the basis of the predicted evolution of the mass inflow rate, this transi- 
tion occurs at At ~ 15 - 21 d. The At = 40 d SED can be modelled 
as arising from a cool ring (T ~ 4500 K) with a large inner radius. 
The At = 180 d SED can be interpreted as the emission coming from 
a larger and colder ring (T ~ 650 K, resembling a protoplanetary 
disk) plus a contribution from the neutron star heated by the accreted 
material (Supplementary Information). 

Estimating the rate of such events is challenging. A wandering 
neutron star can intercept minor bodies as it passes through a plan- 
etary system. The capture rate can be evaluated on the basis of the 
hypothesis that all stars have Oort clouds similar to that of the Sun 25 . 
Applying this to the case of GRB 101225 A, we derive a rate of ~ 0.3 
events per year, in line with Swift observations. However, indirect ar- 
guments suggest a lower capture rate 26 and, in addition, larger minor 
bodies are rarer. The likelihood of a neutron star retaining its origi- 
nal population of small bodies is rather low, because they are unlikely 
to survive the supernova event. Nevertheless, planetary-like systems 
can reform around millisecond radio pulsars 27,28 , providing a viable 
mechanism for the occurrence of GRB 101225 A. 
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Figure 1 | Light curves of GRB 101225A. GRB 101225A light curve in five 
energy bands: X-rays at 1 keV (black dots), UV at 2,030 A (green) and 
2,634 A (blue), and optical at 6,400 A [R band, red) and 7,700 A (I band, 
orange). Error bars, 1 a. The black arrow indicated the XMM-Newton 
upper limit, other arrows optical upper limits. The X-ray light curve only 
refers to the disk contribution to the total flux (~ 0.3 of the total, as de- 
rived from spectral modelling) and is corrected for the interstellar absorp- 
tion column density, of N H = (2.0 ± 0.1) x 10 21 cm~ 2 (which is greater 
that the Galactic value, 7V° al = 7.9 x 10 20 cm -2 ). This enhanced col- 
umn density may be due to the fraction of the minor body mass that has 
been expelled from the system. The continuous lines of different colors 
(the same as the data) represent the fit to the light curves using the phe- 
nomenological model of tidal disruption 21 . Because the model predicts 
a late transition of the accretion disk to a 'cold' solution, the fit has been 
carried out up to ~ 20 d (the excluded region is indicated in pink-gray). 
A thick pink line indicates the time of the transition. Our model has four 
parameters: the mass of the minor body (A/», we assume for simplicity a 
density of 1 g cm~ 3 , thereby fixing the radius, i?»), the periastron (r p ), a 
geometrical factor (D 2 / cosi, where D is the source distance and i the 
source inclination) and the optical absorption (A v ). The best fit is ob- 
tained for: M, ~ 5 x 10 2n g, r p ~ 9 x 10 3 km, D/Vcosi ~ 3 kpc 
and A v ~ 0.8 (in excess of the Galactic A° al = 0.3, consistent with 
the value determined by X-ray analysis). The peak mass accretion rate 
with these parameters is M « 2 x 10 16 g s _1 and the peak luminosity 
L ?s 3 x 10 36 erg s~ 1 , consistent with our hypothesis of a sub-Eddington 
accretion. In this regime, no emission lines are expected, at variance with 
super-Eddington events. Inset: X-ray light curve as observed by the Swift 
X-Ray Telescope. Shaded regions highlight the periastron passages cal- 
culated using our tidal disruption model, which have a timescale of 6,000- 
10,000 s. Outside each band, the X-ray emission markedly decreases, 
and this effect is strongest for the first passage. We verified that this effect 
is not instrumental in nature nor due to spectral changes. Inhomogeneities 
in the return debris are expected to give rise to variability on the fall-back 
timescale. 
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A Christmas minor body falling onto a neutron star: 
Supplementary Information 

Supplementary Methods 

X-ray data analysis 

Swift XRT & UVOT data analysis The X-ray Telescope (XRT 28 ) began observing GRB 101225A on 2010 De- 
cember 25 at 9:00:48.5 UT, 1383 seconds after the BAT trigger 29 , and ended on 2011 January 04 at 12:20:38 UT, 
collecting 3 ks in Windowed Timing (WT) mode, and 120 ks in Photon Counting (PC) mode, spread over 11 days. 
Supplementary Table 1 reports the log of the XRT observations. 

The GRB 101225A light curve at 1 keV shown in Fig. 1 of the main text has been obtained with the Swift burst- 
analyser provided by the UK Swift Science Data Centre at the University of Leicester 30 . The original curve, already 
corrected for absorption, has been multiplied by a factor of ~ 0.3 to obtain the flux at 1 keV from the flux in the 0.3-10 
keV energy band, obtained from spectral modeling (see below). 

We extracted the XRT spectra of the first five orbits covered by Swift for spectral energy distribution (SED) pur- 
poses. The XRT data were processed with standard procedures (XRTPIPELINE vO.12.6), filtering and screening 
criteria by using FTOOLS in the HEASOFT package (v.6.10). Both WT (first two orbits) and PC events were consid- 
ered. The selection of event grades was 0-2 and 0-12 for WT and PC data, respectively 28 . We corrected for pile -up 
when required. The spectra were also corrected for point spread function losses and vignetting generating the ancil- 
lary response files with XRTMKARF. We used the latest spectral redistribution matrices (v.012). The XRT analysis was 
performed in the 0.3-10 keV energy band. Data from the first two orbits were taken in WT mode and spectra were 
rebinned to have at least 50 counts per energy bin. Data from the last three orbits were obtained in PC mode and were 
rebinned in order to have 20 counts per energy bin. 

The Ultraviolet/Optical Telescope (UVOT 31 ) observed the target simultaneously with the XRT. A first account 
of these observations can be found in ref. 32. We extracted two reference light curves in the UVW1 and UVW2 
bands. The data analysis was performed using the UVOTIMSUM and UVOTSOURCE tasks included in the HEASOFT. 
The latter task calculates the magnitude through aperture photometry within a circular region and applies specific 
corrections due to the detector characteristics. The reported magnitudes are in the UVOT photometric system 33 . 

We also extracted five SEDs coincident with the first five orbits of the Swift satellite. For the first two orbits data 
in all filters were available (UVW2, UWM2, UVMl, U, B and V), while for the remaining three only V, B and UVW2 
were present. We summed different exposures in each orbit, and extracted a spectrum with the task UVOT2PHA for 



Supplementary Table 1 | Swift XRT observing log. 


Obs.ID. 


Start time 


End time 


XRT exp. time (s) 


00441015000 


2010-12-25 UT18:21:55 


2010-12-26 UT03:01:15 


12527 


00441015001 


2010-12-26 UT02:57:31 


2010-12-27 UT23:47:25 


33093 


00441015003 


2010-12-28 UT01:06:00 


2010-12-28 UT16:37:09 


9694 


00441015004 


2010-12-29 UT01:08:01 


2010-12-29 UT16:41:51 


10895 


00441015005 


2010-12-30 UT01:08:01 


2010-12-30 UT16:40:26 


11033 


00441015006 


2010-12-31 UT02:49:01 


2010-12-31 UT23:00:15 


9576 


00441015007 


2011-01-01 UT01:27:01 


2011-01-02 UT00:41:47 


9949 


00441015008 


2011-01-02 UT01:28:01 


2011-01-02 UT23:25:17 


11137 


00441015009 


2011-01-03 UT00:04:32 


2011-01-03 UT23:30: 15 


10853 


00441015010 


201 1-01-04 UT06:39:01 


201 1-01-04 UT12:20:38 


4172 
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Supplementary Table 2 | Swift XRT and UVOT SEP model fitting. 



Model* 


Red y 2 (Actf\ 


N.h.p. 


+ 


BB+DiskBB 


1.15 (848) 


2x 10 


-3 


BB+PL 


1.69 (848) 


2x 10~ 


-32 


BB+Brems. 


2.82 (848) 


1 x 10" 


146 


DiskBB+PL 


1.36 (848) 


2x 10" 


-11 


NSA+DiskBB 


1.15 (852) 


2x 10 


-3 


NSA+PL 


1.45 (852) 


2x 10 


16 


DiskBB+COMPTT § 


1.04(842) 


0.17 





* Models adopted in the fit of the five SEDs described in the text. The model consists of two components (BB= black 
body; DiskBB= emission from an accretion disk; PL= Power Law; Brems.= bremsstrahlung emission; NSA= Neutron 
Star Atmosphere with nsatmos spectral model). All the double components model include absorption (using the 
model TBABS). Spectral fits were carried out using the XSPEC package. 
+ Null hypothesis probability. 

§ The normalisation N of the disk black body component is relatively small. The relation between the disk inner radius 
r; n and N should be r m x V cos i = y/N x Dio km (where i is the disk inclination and Dio the source distance in 10 
kpc units). However, the correct prescription to transform the normalization into the inner disk radius depends on the 
effective temperature to color temperature ratio (which is basically unknown in our case 34 ). In any case the relatively 
low value of the normalization calls for a relatively high disk inclination. 

each filter and orbit. The extraction region was chosen to be 6 arcsec radius for UV filters and 3 arcsec radius for optical 
filters due to a contaminating source. A separated closeby background region with 14 arcsec radius was adopted. 

Early time spectral energy distribution We fit together the five data sets within the spectral energy package 
XSPEC. The spectra cannot be fit with any single (absorbed) component model. We thus consider a two-component 
model (see Supplementary Table 2). For all the SEDs we keep tied the absorbing column density (modeled with 
TBABS 35 ) to the same value. The Galactic column density is A^ al = 7.9 x 10 20 crrr 2 (ref. 36) and the galactic 
extinction is A^ al <~ 0.3 (ref. 37). A model made by a black body (BBODYRAD) emission plus the emission from a 
multicolor accretion disk (DISKBB) provides the best fit to the data. Alternatively, also the neutron star atmosphere 
(NSATMOS 38 ) model plus a disk black body emission provide a comparable good fit. 

These two models can be easily interpreted in the present framework as emission from the accretion disk and from 
the neutron star surface heated by the accreting matter. The black body model and the neutron star atmosphere models 
can be regarded as two extremes of the emission from the neutron star surface. The best fit model is made by a disk 
component plus a black body component (or atmosphere component). The presence of this latter component requires 
the presence of a hard surface. This clearly rules out black holes. 

In the framework of accretion onto a neutron star surface, the black body or the neutron star atmosphere mod- 
els represent a first attempt to model the physics of the process. Matter accreting from the disk needs to reach the 
more slowly rotating neutron star. This happens in the so-called boundary layer where a large fraction of the accre- 
tion energy is released. The boundary layer region is quite complex. The gas can reach very high temperatures, so 
Comptonization of soft photons is (likely) the dominant energy loss mechanism 39,40 . Boundary layer spectra are fit 
with Comptonization models, with COMPTT (within XSPEC) being the usual choice 41 . We fit the same SEDs with 
a model made by a multicolor accretion disk (DISKBB) plus a Comptonization model (COMPTT), and an absorption 
component (TBABS). The D I SKBB normalization is tied together in all the SEDs, depending on the disk inner radius. 
The fit is definitely better than the one with the black body plus accretion disk, with a reduced = 1 .04 (842 degrees 
of freedom, d.o.f.), although the disk parameters as well as the overall luminosities are similar. Spectral parameters 
are reported in Supplementary Table and in Supplementary Figure 1. The source unabsorbed luminosity in the 0.3-10 
keV energy range is 1.7 x 10 36 , 4.3 x 10 35 , 1.9 x 10 35 , 1.1 x 10 34 and5.9 x 10 34 erg s _1 , respectively. The bolometric 
corrections estimated over the 0.0001-1000 keV interval are 2.2, 1.8, 2.3, 2.3 and 1.6, respectively. 
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Supplementary Table 3 | Spectral parameters of the Componization plus multicolor accretion disk models of the 
five SEDs. 



Mid. time 


Disk temperature 


Seed Temperature 


Plasma Temperature 


Optical depth 


(hr) 


(keV) 


(eV) 


(keV) 




0.43 
1.70 
3.34 
4.95 
6.55 


1 C1+ 016 

1 - 81 -0.18 

1 oo +0 - 10 

1,uu -0.08 

<0.59 
<0.53 
<0.48 


49+^ 
17+1 
14+1 
10+J 
lit] 


86+- 
56+} 

61^ 

J ->-41 

10+J 


0.38+1!;- 
0.50+°;}? 
51+ - 20 
63+ 30 

u - o:, -0.11 

<2.08 



The column density isN H = (2.0 ±0.1) x 10 21 cm" 2 . 

All the errors are 90% confidence level for one parameters of interest, i.e. have been computed for a Aj 2 = 2.71. 




Supplementary Figure 1 | Evolution of the spectral energy distribution with time. Swift XRT and 

UVOT unfolded spectra taken during the first five orbits were used to build spectral energy distributions 
(SED). The mid-times of the SEDs are 0.43 (black), 1 .70 (red), 3.34 (green), 4.95 (blue), 6.55 (magenta) 
hr, respectively. Here we show these five SEDs based on the best fit model made by a multicolor disk 
black body and a boundary layer emission model. Their contribution is split and indicated with dotted 
lines. This model provides a fair description of all the data with a reduced x, 2 ed = 1.04 (842 degrees of 
freedom, d.o.f.). On the contrary, a black body plus a power law component model provides a reduced 
Xred = 1-69 (848 d.o.f) and a disk black body plus a power law component model provides a reduced 
Xmd = 1-36 (848 d.o.f). These are firmly ruled out. Given the upper limit derived for the disk temperature 
in the last three SEDs, the disk component shown here is just for the best fit temperature and does not 
reflect the maximum disk component that it is possible to include in the data. In the first two spectra 
the multicolor disk black body component comprises ~ 30% and — 10% of the total flux. The source 
unabsorbed luminosity in the 0.3-10 keV energy range is 1.7 x 10 36 , 4.3 x 10 35 , 1.9 x 10 35 , 1.1 x 10 34 
and 5.9 x 10 34 erg s -1 , respectively. 
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Sharp drop of the X-ray flux One of the key predictions of our model is that accretion of matter in the early stages 
of the minor body disruption is not steady and non-homogeneity in the retuning debris are expected during the first 
periastron passages. These give rise to variability on a fall-back timescale, providing a quasi-periodicity (see e.g. Fig. 
10 in ref. [22]). 

Signs of this variability are easily visible in the X-ray light curve (see inset in Fig. 1). We focus here on the 
data in the second orbit in WT mode (4940-7270 s). We positively check that the satellite attitude is stable and the 
source is far from bad columns (whose presence is anyway accounted for and corrected). We select the first 165 s 
(1272 counts) during the flux rise and the following 484 s (6286 counts). We build the appropriate exposure maps 
and ancillary response files and fit the two spectra with the same absorbed power law with only the normalization free 
to vary. The fit provides an acceptable description of the data with a reduced # r 2 d = 1.00 (with 274 d.o.f. and null 
hypothesis probability of 50%). This indicates that the variation is not induced by a spectral change, either an increase 
of the absorbing column density or a change in the power law photon index. This change can instead be interpreted as 
variation in the mass accretion rate producing a flux variation of the same amount. 

We repeat the same check on the late curve. We took 276 s (3325 counts) and 103 s (704 counts) before and after 
the flux decay. We repeat the same fit with an absorbed power law, finding a reduced j r 2 ed = 1.06 (with 159 d.o.f. 
and null hypothesis probability of 27%). Also in this case the variation is consistent with the hypothesis of not being 
induced by a spectral change. 

XMM-Newton data analysis XMM-Newton 42 observed the field of GRB 101225A starting on Jan. 17,2011 UT18:16:32 
(23 d after the burst) for <~ 32 ks under a Directory discretionary Time program (PI Campana). All the three EPIC 
cameras were operated in Full Window mode with thin filters. We generated final product using SAS v. 9.0.0, using 
the latest calibration files. A high background interval affected the end of the observation resulting in 23 ks of effective 
exposure for MOS1 and MOS2 cameras, and 21 ks for the pn camera. No source is detected at the position of GRB 
101225A in all images. Assuming a worse-case power law spectrum with photon index F = 2 and the Galactic column 
density, we derive a 3(7 upper limit on the 0.5-10 keV unabsorbed flux of <~ 10~ 14 erg cm~ 2 s -1 . A black body 
spectrum with temperature 0.3 keV, would result in a factor of ~ 2 tighter upper limit. 

Optical/UV data analysis 

TNG data analysis We imaged the field of GRB 10 1225 A with the Italian 3.6m Telescopio Nazionale Galileo 
(TNG), located in La Palma, Canary Islands. Optical R and /-band observations were carried out with the Device 
Optimized for the LOw-Resolution (DOLoReS) camera on 2010 Dec. 31, 2011 Jan. 16-17, Feb. 7 and Feb. 17-19. 
A /-band NIR observation was performed with the Near Infrared Camera Spectrometer (NICS) on 2011 Feb. 8. All 
nights were clear, with seeing in the range 0.9" — 1 .2". The complete observing log is reported in Supplementary Table 
4. Image reduction was carried out following the standard procedures: subtraction of an averaged bias frame, division 
by a normalized flat frame. NIR frames were reduced using the jitter pipeline data reduction, part of the ECLIPSE 
package A . Astrometry of the optical and NIR images was performed using the USNOB1.0 B and 2MASS C catalo gues, 
respectively. Aperture photometry was made with the ESO-MIDAS D DAOPHOT task for all objects in the field. The 
photometric calibration was done against Landolt standard stars for optical images and against the 2MASS catalogue 
for NIR images. In order to minimize any systematic effect, we performed differential photometry with respect to a 
selection of local, isolated, and not-saturated reference stars visible in the field of view. 

Optical light curve To build the light curve shown in the main text several observations were taken from Gamma- 
ray Burst Coordination Network (GCN) circulars. In Supplementary Table we give a full account of them. Other 
observations were taken from Thone et al. 11 , reporting GCT observations. 

A http://www.eso.org/projects/aot/eclipse 
B http://www.nofs. navy.mil/data/fchpix 
c http://www.ipach.caltech.edu/2mass 
D http://www.eso.org/projects/esomidas 
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Supplementary Table 4 | TNG Optical/NIR observing log. 



Time of obs. 


A T 


Telescope 


Filter 


Exposure 


Magnitude 


Flux 


(UT) 


(d) 






(s) 


(Vega) 


(MJy) 


20101231.86502 


06.08880 


TNG (DOLoReS) 


R 


3x300 


23.75±0.13 


0.91 ±0.11 


20110116.84203 


22.06581 


TNG (DOLoReS) 


R 


7x300 


24. 10 ±0.29 


0.66±0.18 


20110117.85127 


23.07505 


TNG (DOLoReS) 


R 


12 x 300 


24.74 ±0.43 


0.37±0.15 


20101231.87893 


06.10271 


TNG (DOLoReS) 


I 


3x300 


23.36 ±0.26 


1.03 ±0.24 


20110207.87314 


44.09692 


TNG (DOLoReS) 


I 


20 x 300 


23.64 ±0.24 


0.80±0.17 


20110217.86948 


54.09326 


TNG (DOLoReS) 


I 


15 x 120 






20110218.85349 


55.07727 


TNG (DOLoReS) 


I 


20 x 120 






20110219.84790 


56.07168 


TNG (DOLoReS) 


I 


13 x 180 






20110218.85696* 


55.08074 


TNG (DOLoReS) 


I 


6540 


>24.5 (3(7) 


<0.36 (3 a) 


20110208.86950 


45.09328 


TNG (NICS) 


J 


56 x 15 


>20.7 (3(7) 


<5.37 (3(7) 



* Average of Feb. 17, 18 and 19 observation epochs. 
Errors at the 1 a level. 



Supplementary Table 5 | Optical observations from GCNs. 



Time of obs. 


Magnitude 


Refs. 


Time of obs. 


Magnitude 


Refs. 


(hr after T Q ) 


(Vega) 


(hr after 7b) 


(Vega) 






fl-band 






/-band 






1.14 


> 21.5 


43 








1.59 


22.1 ±0.1 


44 


1.59 


21 .5 ±0.2 


44 


5.47 


22.6 ±0.2 


45 


5.47 


21.6±0.2 


45 


7.08 


22.6 ±0.2 


46 


7.07 


22.3 ±0.2 


13 


7.22 


22.2 ±0.2 


13 








24.2 


22.2 ±0.2 


47 


24.2 


21.6±0.2 


47 


28.2 


22.8 ±0.2 


48 


28.5 


21 .7 ±0.3 


13 


50.0 


22.6 ±0.2 


13 








55.2 


22.8 ±0.2 


49 


55.2 


22.3 ±0.2 


49 


123 


23.3 ±0.2 


50 








242 


24.1 ±0.1 


51 


244 


23.6±0.2 


13 


338 


24.5 ±0.2 


52 








432 


24.4 ±0.1 


14 








512 


24.2 ±0.1 


13 








674 


24.6 ±0.1 


13 








939 


25.1 ±0.1 


13 


952.3 


24.3 ±0.2 


13 


972 


25.0 ±0.1 


13 


952.4 


24.2 ±0.2 


13 


1008 


> 26.0 


53 


1486 


24.8 ±0.3 


13 



Errors at the 1 a level. 
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Modelling of the light curves: tidal disruption model 

We use the tidal disruption model of Lodato & Rossi 21 , scaled down to a system comprised of a minor body object 
being disrupted by a neutron star. We also restrict the model to sub-Eddington luminosities, where the wind component 
- considered in [21] - is absent. 

We thus consider a minor body of mass M % and radius /i,, on a parabolic orbit around a neutron star of mass 
Mns = 1 -4M Q and radius R^s = 12 km, with periastron equal to r p . For simplicity, we assume that the minor body has 
a density p = 1 g cm~ 3 . If the periastron is smaller than the tidal radius r t <~ 10 5 — 10 6 km, the minor body gets tidally 
disrupted. The presence of a magnetic field (if any) does not affect the infalling body until matter is fully ionized. 
Nearly half of the debris become unbound, while the other half are launched into highly eccentric orbits. We follow 
here the Rees' 1 suggestion, that assumes a flat distribution of mechanical energy in the debris. Taking into account 
the internal structure of the object may lead to some correction to this model at early times 54 . The most bound debris 
return to periastron after a time f m j n , given by: 



3/2 / r 3 



(1) 



min 21/2 y R j y GMns 
The other debris falls on a longer time depending on the specific energy E, with a fall-back time of 

2nGM NS 

ffb ~(-2£)3/2 (2) 

Clearly the fall-back time has a strong dependence on the periastron radius (to the third power) and small changes in 
r p , determine large changes in f m j n and to,. Based on our model, we can estimate that f m j n is in the 6,000-10,000 s 
range. 

In the following, we identify f m ; n with the time of the BAT trigger. After f m ; n , the rest of the debris return to 
periastron at a rate: 

J 'nun \ f min / 

The debris shock at periastron and rapidly circularize to form a thin accretion disk. The viscous timescale in the 
disk is initially shorter than the fall-back timescale and the disk can thus be approximated by a sequence of steady 
state models, extending from R{ n ~ Rns to the circularization radius R om = 2r p . The effective temperature of the disk 
as a function of disk radius and time is given by 

^( S ,,)=-M[.-,/^|. (4) 

The emitted flux at frequency v is computed assuming a multi-colour blackbody disk spectrum: 

cos; f R ^2hv 4 IkRAR 
VFv[t) --&J Rm C 2 ehv/kmt) _ V ( 5 ) 

where D is the distance to the source and ;' is the inclination of the disk with respect to the line of sight. We then redden 
the optical and UV fluxes, assuming the Cardelli 55 reddening law, with a visual extinction coefficient Ay, taken as a 
free parameter. 

In total, our model has thus four free parameters: the mass of the minor body M», the periastron r p , the geometrical 
factor D 2 / cos; and A v . We find that the following set of parameters provides a good match to the X-ray, UV and 
optical light curves: M* = 5 x 10 20 g (which corresponds to a radius R* w 50 km for the selected density), r p = 9x 10 3 
km, D/Vcos/ = 3 kpc and Ay = 0.75. Note that the resulting periastron is well within the expected tidal radius for a 
minor body, and therefore consistent with the hypothesis that the object is tidally disrupted. The radius of the minor 
body is somewhat large but if the density is larger the radius will shrink. For a density of 10 g cm~ 3 the radius will be 
pa 24 km. Despite this the involved mass is not small in the framework of minor bodies. 
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The peak mass accretion rate with these parameters turns out tobeM(f m ; n ) k2x 10 16 g s -1 and the peak luminosity 
is thus L(? m j n ) ~ 3 x 10 36 erg s -1 , consistent with our hypothesis of a sub-Eddington accretion. Note that this will imply 
that no emission lines are expected from super-Eddington wind as expected in tidal disruption models. Comparing this 
luminosity with the BAT peak flux (including a factor of 2.2 for the bolometric correction based on spectral modeling) 
we obtain a distance of ~ 3.5 kpc. We assumed a distance of 3 kpc, as motivated also by the matching of the observed 
fluence with the conversion of the minor body gravitational energy into radiation. The two estimates match with very 
good approximation. This distance places the neutron star in the Perseus arm. Its height above the Galactic plane is 
~ 1 kpc. 

Late time spectral energy distribution 

It is well known that accretion disks around compact objects are thermal-viscous unstable at a temperature corre- 
sponding to partial ionization of hydrogen. A model based on this instability captures the main properties of dwarf 
nova and transient low mass X-ray binary outbursts 22 . The same physics applies to disk made of helium or of heav- 
ier elements 23 , leading to a critical mass inflow rate at a given radius, where the instability sets in and then quickly 
propagates to the entire disk. In our case, the accretion rates decreases with time (Eq. 3), the disc cools down and 
eventually will hit the critical ionization temperature. The outer disc radius is the first radius at which this occurs. In 
our sub-Eddington regime, the outer disc is always gas pressure dominated, where the main opacity is Kramer (e.g. 
Fig 5.4 in ref. 56). Under these conditions, the temperature at outer disk R out scales with the accretion rate as 

T(R out ) = 9.6 x 10 3 (M_ 13 (f)) 3/10 K, (6) 

where M_i3 is the mass accretion rate in units of 10~ 13 M Q yr -1 and we assumed a standard Shakura & Sunyaev 
disk and a viscosity parameter a = 0.1. For low metallicity disc, the critical accretion rate is around 10~ 13 M yr -1 
(corresponding to transition temperature of T <~ 10 4 K, due partial hydrogen ionization 24 ). From Eqs. 3 and 6, we 
thus get that the the transition happens at t ^ 15 d. It is well-known that the critical temperature decreases with 
metallicity 24 . Our disc will likely have solar metallicity or higher, so 15 d is just a lower limit. If we take a slightly 
lower ionization temperature of <~ 8,000 K, we obtain a transition at t ~ 20 d. We thus conservatively assume a 
transition time around 20 d. After the transition to the 'cold' phase the disk emission would drop faster in the blue 
optical filters and slower in the red optical filters. 

A second known effect occurs to the disk for decreasing mass inflow rates. For high mass inflow rate the mass 
inflow can reach the neutron star surface developing a boundary layer where matter is slowed down before accreting 
onto the neutron star surface. If the neutron star possesses a significant magnetic field (10 8 — 10 9 G), for a sufficiently 
low mass inflow rate this magnetic field will be able to disrupt the disk flow at a magnetospheric boundary. This 
happens when the magnetic pressure CP mag (r) °< B 4 r~ 6 , with B the neutron star magnetic field) equates the disk ram 
pressure (PdiskM Mr -5 / 2 , where M is the mass accretion rate) at a radius larger than the neutron star radius 56 . 
When this occurs the disk inner edge ends at the magnetospheric boundary and for smaller radii the motion of the 
infalling matter is controlled by the neutron star magnetic field. Since the magnetospheric radius is proportional to 
rM x B 4 / 1 ' M~ 2 I~' ', larger magnetospheric radii are expected for lower mass inflow rates. Since M decrease as f~ 5 / 3 , 
being a tidal disruption event, we can impose to our model that the inner radius of the disk r m {t) t~ w l 2X . 

With these well-know ingredients we approached the late time SEDs. At very late times ( £ 20 d), the mass accre- 
tion rate is very low and we expect a large magnetospheric radius (if the neutron star has a non-negligible magnetic 
field). Given that the outer disk is ~ 2 r p ~ 2 x 10 9 cm, we approximate the disk in this time interval as a small ring and 
treat its emission like a black body. We fit an optical-IR SED at 40 d taken from ref. 1 1 . A black body fit (absorbed at 
Ay ~ 0.75 as derived from the light curve fitting) provides a good description of the data, with a reduced % 2 = 1 .7 for 3 
d.o.f. (16% null hypothesis probability). The resulting temperature and radius are 7bb ~ 4500 K and /?bb ~ 2.2 x 10 8 
cm (at 3 kpc, see Supplementary Figure 2). The derived temperature and emitting radius nicely fit with the predictions 
of our late time disk model. 

Taking /?bb as the radius of the magnetosphere at 40 d, based on the model predictions about the evolution of the 
mass accretion rate we can estimate the magnetic field of the neutron star, which is B ~ 10 9 G. This is in line with the 
magnetic field of millisecond radio pulsars. 

A very late time detection (At = 180 d) has been reported 11 . These data show a large decrease in the r 1 band from 
data taken at Af = 40 d (factor of 4.6 in flux) but a mild decrease in the bluer g' band (factor of ~ 1.5 in flux). To 
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Supplementary Figure 2 | Late time (A/ = 40 d) spectral energy distribution. The late time spectral 
energy distribution taken from Thone et al. 11 has been fit with a absorbed (Ay = 0.75 as derived from 
the light curve fitting) black body emission model. The best fit model provides a temperature of 4500 K 
and a radius of 2.2 x 10 8 cm (at 3 kpc). The reduced j r 2 ed = 1.7 with 3 degrees of freedom, implies a null 
hypothesis probability of 16%. 



model this spectrum we need a hot component that we identify with surface emission from the neutron star heated 
by the accretion episode and a cold component that may come from a dusty ring/disk, where the minor body formed. 
Assuming a neutron star with a temperature of ~ 10 6 K we obtain a disk inner radius of ~ 10 13 cm and a temperature 
of ~ 650 K. Clearly with just two points this is just indicative E . 

Why it cannot be another object? 

"When all possibilities have been eliminated, whatever remains, however improbable, must be the truth" (Sir Arthur 
Conan Doyle). 

Long-duration Gamma-ray Burst GRB 101225A has been discovered as an image trigger. Inspection of the BAT 
light curve indicates that there has not been a clear peak but a fiat emission 10 . The event has also been detected by 
MAXI/GSC in the 2-10 keV energy band 57 . The burst is extremely long and the X-ray spectrum collected by the BAT 
is soft. These peculiar characteristics resemble those of GRB 060218 13 and GRB 100316D 58 . Also the UVOT flux 
is bright and similar to the shock break out emission seen in GRB 060218, but not in GRB 100316D. Moreover, at 
variance with GRB 060218 and GRB 103016D, GRB 101225A does not show signs of a supernova (SN). This might 
not be an insurmountable problem since other two low energetic and close by (z < 0.13) GRBs without SN have been 
observed in the past 59-61 . One of them, GRB 060614, has also a relative long duration (790 = 102 s, ref. 62), whereas 

E Assuming the same neutron star component to be present in the At = 40 d data, the best fit parameters for the disk change only mildly to 
#BB ~ 2.8 x 10 s cm and 7bb ~ 4000 K 



8 



the other (GRB 060505) might be a short GRB. GRB 060614 had a quite spiky and hard X-ray emission (peak energy 
£p > 100 keV) as detected by the BAT and it had a bright X-ray afterglow. In addition, GRB 060614 had a fast UV 
emission too but its emission declined achromatically 63 . 

All these bursts are very bright in XRT allowing us to disentangle the absorption pattern from our Galaxy (of a 
fixed amount) from that within the host galaxy, providing the redshift of the GRB 64 (see also S. Campana et al. 2012, 
in preparation). The redshift of GRB 1003 16D has been predicted to be 0.014 < z < 0.28 (90% confidence level 65 ), in 
agreement with the observed value of z = 0.059 (ref. 58). We applied the same technique to GRB 060614 providing 
0.12 < z < 0.20 (measured z = 0.125) and to GRB 060218 providing 0.01 < z < 0.08 (measured z = 0.033). 

Fitting the XRT X-ray spectrum of GRB 101225 A with a Galactic plus intrinsic host galaxy absorption power law 
plus black body model led to the determination of the redshift z = 0.07+q ^ (90% confidence level) with a reduced 
X 2 = 1.07 (505 d.o.f. 66 ). Adding to these data also the optical/UV data described above and keeping fixed the Galactic 
column density to the measured value, the fit results to be highly unacceptable (x r 2 d = 2.01, 847 d.o.f.). A fit with a 
broken power law (instead of power law) with a difference in the power law indices of 0.5, as sometimes done when 
fitting such a large energy range based on GRB afterglow theory, results in an unacceptable fit too (j r 2 d = 1-45. 842 
d.o.f.). This clearly points against a GRB whose X-ray spectrum is always described a power-law. 

Moreover, a low redshift solution has also to face the problem of a very faint optical galaxy. A host galaxy 
candidate has been found in Gran Telescopio Canarias (GTC) data with g' AB <~ 27.2 and r' AB ~ 26.9 (ref. 11). At a 
redshift of z = 0.07 any host galaxy beneath the transient would have to be very faint (M g <~ —9, i.e. fainter than a 
large globular cluster, e.g. co Cen). To make the host galaxy and the SN brighter one could think to move the object 
further away, but still within z 1 . 1 due to the detection in the UVW2 band. In this case however, the energetics of the 
GRB will become much larger and it will be a strong outlier of the Amati relation 67 . To lie within 3 a from the Amati 
relation best fit, the maximum allowed redshift is z ~ 0.3. This would change the distance modulus by just <~ 3.5 mag, 
still leaving an extremely faint host galaxy (see also below). 

Given all these difficulties, we regard the GRB nature for GRB 101225A as highly unlikely. 

Short-duration Gamma-ray Burst A short duration GRB might mimic GRB 101225A if we are observing the so- 
called extended emission phase 62 . The problem with this interpretation is that no hard X-ray spike has been observed 
associated with this event. Even if for some reason this spike has been overlooked, the extended emission observed in 
several short GRBs is too short compared (lasting only 100 to 1,000 s, ref. 68) with what observed in GRB 101225 A, 
and no bright UV long-duration emission has ever been detected in short GRBs. We therefore exclude this possibility. 

Bright Supernova Several cases of association between core-collapse SNe and GRBs have been discovered 69 . 
Given the large luminosity at maximum of these SNe, the association between GRB 101225A and a bright SN is 
obviously ruled out by the lack of bumps in the optical light curve. Assuming a SN template light curve and moving 
it at several redshifts (including bolometric correction) we can exclude that a SN1998bw-like and a SN2006aj-like are 
present in the data for z 0.8. 

Faint Supernova A different possibility is that extremely sub-luminous GRBs may be associated with ultra-faint 
Core-Collapse Supernovae 70 ' 71 . Indeed, Thone et al. 11 suggested this possibility to explain the peculiarities of GRB 
101225A. They fit the late time SED with a SN1998bw template spectrum finding a good matching with a <~ 10 fainter 
supernova. This result has been obtained by varying contemporaneously the SN redshift (found to be z = 0.33t_Q Q 1 4 ) 
and the stretching factor (which is determined for all the other SNe fixing the redshift to the known value). We also 
note that all the other templates used (except two) provide statistically acceptable fits, spanning the redshift range 
0.2-0.5. 

We fit the late light curve (> 100 hr. i.e. ~ 4 d) in the R and / bands (excluding the very late detection at ~ 180 
d). While in both light curves a small bump could be identified, a fit with a simple power law provides a statistically 
sound description of the data. In the R band we obtain a reduced # r 2 d = 1 .0 (9 d.o.f., see Supplementary Figure 3) and 
in the / band we obtain a reduced j r 2 d = 0.9 (5 d.o.f.). Thus, the presence of a bump is not justified from a statistical 
point of view. Note also that the bump, if present at all, should peak at different times in the two bands: at 22 ± 6 d 
(90% confidence level) in the R band and ~ 48 d in the / band, respectively. This is due to the fact that in the / band 
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around ~ 20 d no data are available, and in the R band no data are available around ~ 40 d. This fact also testifies that 
the putative bump(s) are just statistical fluctuations. 

A host detection has been claimed at g' = 27.2 ± 0.3 and r' = 26.9 ± 0. 14 (ref. 11). At a redshift of z = 0.33 the 
absolute magnitude is M g = —13.7, that is 2 mag fainter than the faintest host galaxy associated to a GRB-SN (GRB 
060218). This host galaxy is much dimmer than the Magellanic Clouds and has the same absolute magnitude of the 
Fornax galaxy. It is peculiar that the dimmest supernova is also associated to the dimmest host galaxy. 
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Supplementary Figure 3 | Late time R and / band light curves. R (red dots) and / (orange dots) band 
curves after 100 hr from the trigger. The dashed red line represents the best fit power law with index 
-0.8 to the R band data; the dashed orange line the best fit (index -0.4) to the / band data. The R band 
slope is much steeper than the / band, indicating a chromatic behavior. 



Blazar/AGN Blazars show pronounced variability at all wavelengths. However, this variability is also accompanied 
by radio emission, that is lacking in this case. Radio emission has not been revealed with 3 O" upper limit of 60 fiJy at 
5 GHz, 4 d after the trigger 72 , and 42 jiiJy and 30 /iJy at 4.5 and 7.9 GHz, 12 d after the trigger 73 . This excludes the 
blazar option. 

Tidal disruption by Intermediate mass black hole The model used here 21 is self-similar and can provide the same 
results scaling Ml — f 2 M*, M' NS = fM^s, D' = fD, and r' p = / 18 r p . This in principle can allow for an event occurring 
onto a supermassive black hole (requiring a star falling on it) or an intermediate mass black hole (IMBH, ~ 1 , 000 M ) 
in the local group with an object like the Moon (~ 10 25 g) falling on it. A massive black hole is associated to a bright 
galaxy (based on the bulge luminosity - black hole mass correlation 74 ) that is ruled out by the very faint quiescent 
counterpart. Moreover, all the observed events are characterized by much longer timescales (years 75 76 ), including the 
recent case of Swift J164449.3+373451 6 20 . 

The location of GRB 101225A lies somewhat close to the Andromeda galaxy (M31). If GRB 101225A belongs 
to this galaxy it would be in its halo at a distance of at least ^115 kpc. In the distant local group it could also be 
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associated to the dwarf spheroidal galaxy Andromeda XVIII (~ 1.4 Mpc) at <~ 12 kpc from its center . The IMBH 
possibility can be excluded based on the probability of having an IMBH well outside M3 1 or in the outskirts of 
Andromeda XVIII, with a Moon falling on it. 

If the object is indeed extragalactic, the emitted energy is very large. Assuming a redshift z = 0.3 the total emitted 
energy at high energies is > 10 51 (ref. 1 1). This is very large and likely further excludes and extragalactic tidal event. 

X-ray binary transient Matter accreting onto a compact object in a binary system might be another option to 
explain GRB 101225A. In this framework the luminosity is limited to be below the Eddington rate. Assuming a black 
hole with mass 10 M© we derive a distance of D ;$ 75 kpc. This implies that the binary must be of Galactic origin. 
There is a large number of transient X-ray binaries. At 75 kpc, the late time detection from GTC would result in an 
absolute magnitude Mr ~ 7.5. This rules out all giants and supergiants companions (the absorption measured from 
the X-ray column density is low so it will not alter this conclusion). The putative companion should be cooler than 
spectral type K2V. This clearly points to a Low Mass X-ray Binary (LMXRB) system. Outbursts of classical LMXRB 
transients have much longer timescales 77 . Shorter timescales are attained by relatively low luminosity transients 
(named faint transients 78 ) linked to accreting millisecond X-ray pulsars. Binaries misidentified in the past as GRBs 
involve the so-called 'burst-only' sources (sources discovered thanks to Type I X-ray burst activity but with no or 
very weak persistent emission 79 ). The X-ray spectrum during the outburst is usually characterized by a power law 
component that is not present in this case. 

Lowering the burst peak to 10 37 erg s _1 the maximum distance becomes ~ 5 kpc, and the allowed companion 
should be cooler than spectral type M5V. This will also put strong limits on the orbital period of the system to ensure 
Roche lobe contact (the period should be less than a few hours). In such short orbital period binaries the companion 
star is always brighter than the nominal main sequence star, due to irradiation from the primary. Only ultracompact 
binaries do not show this effect (f or b 1 hr). At 5 kpc the XMM-Newton upper limit implies a limiting luminosity 
of <~ 3 x 10 31 erg s _1 . All neutron star transients have a quiescent luminosity larger than this value 80 except one, the 
ultracompact binary candidate HI 905+00 (ref. 81, 82) out of ~ 20. 

Black hole binaries have X-ray quiescent luminosities sometimes lower than the XMM-Newton limit, but in this 
case the formation of a binary system with such a large mass ratio and small orbital period is not straightforward. 
Moreover, the X-ray spectral decomposition calls for a 'surface' component besides the accretion disk. 

The ratio of the X-ray luminosity to the optical luminosity is ~ 10 4 , which is about two orders of magnitude larger 
than observed in LMXRBs 83 . This is accounted for in our model by the small disk whereas in a LMXRB the disk is 
well developed from the truncation radius down to the compact object. 

We also note that the Swift observations during the first two orbits showed a flux decrease at the end of each 
orbit, possibly suggesting a period of <~ 2,500 — 3,000 s. The presence of these dips is not confirmed in following 
observations. In addition, we carefully analyzed the X-ray spectra during these putative dips and outside them during 
the first orbits. Absorption dips in LMXRBs are due to obscuration in the thickened outer regions of an accretion 
disk 84 . A detailed spectral analysis of spectra (as well as a careful check of the satellite attitude) taken during the 
putative dips and outside did not reveal any change in the absorption column density, ruling out this possibility (see 
above). 

On top of all these difficulties there is the position of GRB 101225A in the Galaxy. All X-ray binaries are 
concentrated toward the Galactic center and bulge. The number of LMXRBs within 30° in latitude from the position 
of GRB 101225A is 5 out of 185 binaries in the 4th Low-Mass X-Ray Binary Catalog 85 . 

Given all these constraints we regard the X-ray binary option too contrived. 

Magnetar Magnetars are young (<~ 10 6 yr) isolated neutron stars with a high magnetic field (typically B £ 10 13 
G) powered by magnetic energy rather than spin-down losses 86 . Magnetars emission is usually accompanied by 
repeated short (< 1 s) high-energy flares that are not observed in GRB 101225A at any level. Pulsations are observed 
during outbursts with spin periods in the 2-12 s range. Bright quiescent X-ray counterparts are always observed 
with Lx XL 10 35 erg s _I due to the cooling of the hot neutron star and non-thermal emission associated with the huge 
magnetic field. None of these properties are shared with GRB 101225A. 
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Swift J195509.6+261406 One possible related system is Swift J195509.6+261406, discovered as GRB 070610 by 
the Swift satellite. From this object X-ray and optical flares were observed, even if no pulsations were detected. Given 
the flaring variability and the peculiar afterglow with respect to GRBs it has been suggested that this source can be 
linked to magnetars 87-89 . A Chandra observation performed ~ 2 yr after the main event failed to detect the source in 
quiescence 80 . The 0.3-10 keV unabsorbed flux limit is ~ 10~ 15 erg s _1 cm~ 2 . Tight upper limits were also derived 
in the optical and IR with H > 23, R > 26.0 and i' > 24.5 (ref. 89,87). Given this observational framework Rea 
et al. 80 questioned the magnetar interpretation, based on the too low quiescent X-ray flux. They concluded that the 
source must lie in the Galactic halo to be minimally consistent with the magnetar scenario. Clearly this option is rather 
unlikely given the paucity of magnetars in our Galaxy. Rea et al. 80 suggested instead an X-ray binary origin for Swift 
J195509. 6+261406, even if the very unusual flaring behavior remains unexplained. 

Here we note that more energetic and shorter events than GRB 101225A can occur in the tidal capture of a minor 
body by a neutron star if the tensile strength of the minor body is strong. In this case one can observe a single, 
bright and short event if the object is accreted by the neutron star directly (as was suggested in the past to explain 
short GRBs 9,90 ' 91 . Alternatively if the minor body is disrupted into several pieces still part of an accretion disk, one 
can observe a behavior similar to GRB 101225 A with repeated flares on top. One possibility is therefore that Swift 
J195509.6+261406 is not a magnetar nor an X-ray binary but a tidal disruption event similar to GRB 101225A. 

Nova Recently the 'Nova-like' system V407 Cyg has been detected by Fermi 92 , so in principle one can expect to 
detect an X-ray flare, associated with Nova-like systems. However, this possibility is easily ruled out by the fact 
that Novae at maximum light are extremely luminous at optical wavelengths. For example, fast evolving Novae, on 
time scale of 1-2 d, have an absolute magnitude at maximum of <~ —9 (ref. 93), that corresponds to V <~ 11 — 15 
for novae in the local group of galaxies. These bright counterparts have not been observed in temporal and spatial 
coincidence with GRB 101225 A. Even if one assumes to have missed the maximum, after 15 d a nova in the local 
group of galaxies should still have an apparent magnitude brighter than V ~ 18 (ref. 94). Also this case is excluded by 
optical observations. On average, recurrent Novae are fainter at maximum by 2-3 mag than classical Novae 95 . Current 
optical observations rule out also this possibility. 

Flare star Active stars can give rise to impulsive fast-rising (seconds-minutes) X-ray flares that can last from a few 
minutes up to hours or days. These flares can be orders of magnitude more powerful than those observed in the Sun. 
The brightest flares have been observed in young stars like AB Dor 96 , RS CVn type of stars 97,98 and Algol type of 
stars 99 . Usually these flares are quite soft when compared to a GRB prompt emission, however in a few cases the 
emission has been detected up to 50 keV and beyond (e.g. for Algol 100 , UX Ari 97 , and II Peg 98 ). The energy budget 
observed for the most powerful of these events are in the range 10 36 — 10 37 erg with a peak luminosity of a few 10 32 
erg s _1 . Less powerful (but still more powerful than in the Sun) flares are commonly observed in the UV Ceti-type 
flares stars, with an energy budget of 10 30 — 10 34 erg and a peak luminosity of 10 28 — 10 30 erg s _1 (refs. 101, 102). 
Still they have an X-ray to bolometric luminosity ratio as high as the above type of flare stars. To check if one of these 
type of flare stars could be responsible for the event observed by Swift, we divided the peak luminosity observed for 
the most powerful flares of these stars by the peak flux observed by the Swift XRT and derived the maximum distance 
to which these stars could have been placed (even more stringent values would be derived if we use the Swift BAT 
prompt peak flux). We derived a maximum distance of <, 50 pc for the RS CVn type of stars and < 10 pc of the UV 
Ceti type of stars. At these distances these stars should have an apparent R magnitude of 7-9 for the RS CVn and of 
10-18 (two magnitudes fainter in quiescence) for the UV Ceti. The faintness of the optical counterpart detected during 
the flare and in quiescence by the Swift UVOT and by optical telescopes (including HST) strongly argues against an 
association between a flare star and GRB 101225A. 
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